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ABSTRACT 

We propose a nonlinear self-consistent model of the turbulent non-resonant 
particle acceleration in solar flares. We simulate temporal evolution of the spectra 
of charged particles accelerated by strong long-wavelength MHD turbulence tak- 
ing into account back reaction of the accelerated particles on the turbulence. The 
main finding is that the nonlinear coupling of accelerated particles with MHD 
turbulence results in prominent evolution of the spectra of accelerated particles, 
which can be either soft-hard-soft or soft-hard-harder depending on the particle 
injection efficiency. Such evolution patterns are widely observed in hard X-ray 
and gamma-ray emission from solar flares. 

Subject headings: Sun: flares — acceleration of particles — turbulence — diffusion — 
shock waves — Sun: X-rays, gamma rays 



A solar flare arises due to fast and spatially localized strong energy release and reveals 
itself in electromagnetic radiation and particle flows. Details of the energy release (as well as 
energy storage) are currently debat ed. Models ba sed on the idea of magnetic reconnection are 
the most popular at the moment rtShibatalll999l ). while inter esting alternative possibili ties, 
like balooning instability ( Shibasaki 200ll ). or circuit models ( Zaitsev fc Stepanov 1992 ) are 
discussed as well. 

A common feature of the solar flares (as well as other astrophysical objects with strong 
energy release) is the production of nonth ermal accelerated particles. There are now ample 



evidence of particle acceleration in flares (IMever et al.l 11956 



Cane et al.lll98fil : IChupplll990l : lAkimov et al 



Mathews fc Venkatesan 1990 



19921 ). Accelerated electrons reveal themselves 
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in a variety of nonthermal emissions observed from radio range to gamma-rays. Even rather 
small number of accelerated electrons arising under a weak a ccelera t ion p r ocess can be visi- 



ble wh en the electrons produce coherent radio emission (e.g.. iBenzl ( 119861 ); iFleishman et al. 
( 120021 )). Nonthermal incoherent emissions (gyrosynchrotron and/or bremsstrahlung) are 
detectable when the acceleration is strong enough, providing a considerable fraction of back- 
ground electrons to be accelerated. 

Observations made by Ramaty High Energy Solar Spectroscopic Imager (RHESSI) for 
the recent years have provided us with new st r ingen t const raints on the acceleration mech- 



anism^) operating in flares. iGrigis fe Benzl (120041 . 120051 ) investigated spectral evolution 



during individual subpeaks of the impulsive hard X-ray (HXR) emission and found each 
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concluded accordingly that the ability to reproduce the SHS spectrum evolution of the ac- 
celerated particle population must be an intrinsic property of the acceleration mechanism 
involved; it is the observational property of the acceleration mechanism that is addressed in 
this Letter. 

A number of acceleration mechanisms and mod els have been propo sed to account for 
the particle acceleration in flares (see, for a review, Aschwanden ( 20021 )). Acceleration by 



DC electric fields, both sub-Dreicer and super-Dreicer, has been considered (jHolmanlll985 



Tsunetalll985l ; iHolman fc Benkalll992l ; iLit vinenkol 1 1 9961 ) . This process is able to provide the 
energization of particles up to 100 keV, thus, it can be considered as a possible mechanism 
of pre-acceleration in flares. 



Stochastic acceleration by turbu 



accel e ration in impulsive so l ar flares (IMiller et al 



1992; Petrosian et al. 1994; Park et all 



ent waves is currently assumed to provide the main 
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1997 



1997; iPryadko fc Petrosianl 



Hamilton fc Petrosian 



1998), while the clas- 



sical diffusive sho ck acceleration is believed to play a role in large-scale gradual events. 
Miller et al.l (119971 ) made a detailed comparison between various acceleration scenarios and 
concluded that the stochastic acceleration is intrinsically consistent with the observational 
constraints on the acceleration time, highest particle energy, and the total number of the ac- 
celerated particles. In this Letter we demonstrate that the turbulent stochastic acceleration 
is also naturally consistent with the SHS spectrum evolution of the accelerated particles. 



Grigis fc Bend (20061) noticed that within the standard model of the stochastic accel- 



eration (IMiller et al.lll996l ) the higher level of the turbulence results in harder steady-state 
spectrum of accelerated electrons and vice versa, which looks consistent with the observed 
SHS evolution. It is, nevertheless, unclear if the real evolution of the spectrum will represent 
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the sequence of such steady-state spectra or will behave differently. We note, however, that 
a fraction of stronger events (typically, the proton reach events) displays a different kind of 



the spectral evolution, namely, soft-hard-h arder (SHH) fjFrost fc Dennis! Il97ll ; ICliver et al. 



19861 ; iKipfingerl Il995l; Saldanha et al.ll2008l ). as well as the gradual phase of the impulsive 
events ( Grigis fc Benz 20071 ); we return to this point later. 



Another important point firmly established by RHESSI data analysis (IBrown et al. 



20071 ; iDennis et al.l 120071 ; iHudson fc Vilmerl 120071 ) is that a significant fraction (some tens of 
percent) of the released energy goes i nto nontherma l accel erated particles. This conclusion 
is also confirmed by the radio data. iBastian et al.l (120071 ) performed a calorimetry of the 
accelerated electron energy in a solar flare. They analyzed the radio spectrum evolution of a 
dense flare, when most of the accelerated electron energy was deposited into the coronal loop 
(rather than into the chromosphere). This made it possible to accurately measure the total 
energy deposited by the accelerated electrons into the coronal thermal plasma, which turned 
out to be as high as 30% of the estimated magnetic energy of the flaring loop. These findings 
imply that the back reaction of the acceleration particles on the accelerating agent (e.g., the 
turbulenc e) is not negligible, in full agreement with time-depen dent test particle a nalytical 
solutions (IBykov fc Fleishman! Il992l ) and numerical modeling (ICargill et al.l 120061 ) . so this 
back reaction must be properly taken into account by the acceleration model. 

As we have noticed, the stochastic acceleration of the charged particles is the most 
prom ising candidate for the particle acceleration in flares (e.g., Miller et alll997 ; Grigis Sz Benz 
20061 ). However, the mentioned energetic requirements complicate strongly the whole theory 
of the stochastic acceleration. First of all, the turbulence energy must be large enough to 
supply the accelerated particles with sufficient energy. This me ans that the turbulence is 
strong, so nonlinear effects are important (e.g.. lYan et al.l 120081 ) and the (quasilinear) ap- 
proximation of the weak resonant wave-particle interaction is invalid any longer. Therefore, 
the theory must include a more general transport equation valid in case of strong large-scale 
turbulence. Second, since this turbulence loses a large fraction of its energy to accelerate 
particles, this damping rate must be properly taken into account; thus, one needs to solve 
two coupled equations - one for the particles and the other for the turbulence. 

The renormali z ed theo ry of particle acceler ation by strong turbulenc e was developed by 
Bykov fc Toptygin ( 1990al lbl). see the review by Bykov fc Toptygin ( 1993 ) as well. The form 
of the equation for the averaged distribution function of nonthermal particles depends on 
whether the turbulence is composed of smooth (wave-like) fluctuations o nly or contains 
also the shock fronts and other discontinuities ( Bykov fc Toptygin 19931). Obsery ati 



of the magnetic field complexity in the flare-productive acti ye regions (Abramenkol 12005) 



along with various models of primary energy release in flares (IShibatal ll999 
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2001 



-4- 



Zaitsev fe Stepanovl 11992; IV 



t urbule nce (e.g., 



Miller et al 



1997 



alios 200 7|) suggest many ways of producing long-wave MHD 



Grigis fe Benal2006l ) including possibly a shock wave en- 



semble ( lAnastasiadis fe Vlahod Il99ll ). Although the generation of discontinuities is likely 



under the impulsive energy release (e.g.. lVlahosll2007l ). we limit our consideration to the case 



of smooth long-wave turbulence only, but will include effect of the shocks in a further study. 

We, therefore, adopt the following scenario. The process of flare energy release is ac- 
companied by formation of large scale flows and broad spectra of MHD fluctuations in a 
reasonably tenuous plasma with frozen-in magnetic fields. Vortex electric fields generated 
by the compressible component of the large scale motions of highly conductive plasma will 
result in efficient non-resonant acceleration of charged particles. 

The distribution function N(r,p,t) of nonthermal particles averaged over an ensemble 
of turbulent motions satisfies the kinetic equation 



^_JL X dN _ 1 d ^ 4 dN | 

dt dr a a ^ drg p 2 dp dp % ' 



The particle source term Fi(p) is determined by injection of the electrons and nuclei. Al- 
though we do not consider explicitly the injection process, we note that there are many ways 
to inject particles into the stochastic acc eleration by strong turbulenc e. The possibilities are 
various versions of the DC acceleration ( Litvinenko! 19961 . 2000 . 2003 ) or resonant stochastic 



ac celeration by the small-sca le waves (IMiller et al 



A nice option proposed recently 



by iFletcher fe Hudson! (120081 ) is that the reconfiguration of the preflare magnetic field can 
result in large-scale pulses of the Alfven waves, which in the presence of strong spatial gra- 
dients will generate field-aligned electric regions capable of accelerating electrons from the 
thermal pool up to 10 keV or above. Our further analysis does not depend on the specific 
injection mechanism and details of the particle injection process. Then, the phase space 
diffusion coefficients D and Xa/3 = X $af3 are expressed in terms of the spectral functions that 
describe correlations between large scale turbulent motions (see Bykov & Toptygin, 1993). 
The kinetic coefficients satisfy the following renormalization equations: 



X 



d 3 kdu> 



2T + S 

iu + k 2 x 
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D = - 

9 



d 3 kdu k 4 S(k,u,t) 



(3) 



(2tt) 4 u 2 + k i x 2 ' 

where T(k, u, t) and S(k, uj, t) are the transverse and longitudinal parts of the Fourier compo- 
nents of the turbulent velocity correlation tensor. The equations for T(k, u>, t) and S(k, u, t) 



can be found in iBykovl (120011 ). We reproduce here only the equation for the longitudinal 
spectral function S(k,u,t) responsible for the particle acceleration: 



dS(k,u,t) _ dU s a {k,uj,t) 
dt dk a 



j ST T(k,u,t) - -/ dS S(k,uj,t) - j ap S(k,u;,t). 



(4) 



This full equation includes nonlinear cascading flux, H^{k, u, t), as well as coupling with 
the transverse function (term j ST T(k,u),t)), the true damping (term —j dS S(k,u,t), and 
the damping due to acceleration of the charged particles (term —j ap S(k,u,t)). Although 
all corresponding processes are generally relevant for the turbulence evolution, we found 
that only the last of them is critically important to provide the SHS spectrum evolution 
due to the nonlinear nonresonant particle acceleration by strong turbulence. In case of a 
single scale long-wavelength injection of the turbulent motions (gaussian spectrum with the 
characteristic wave-number ko) we can neglect both cascading term in the left hand side 
and direct turbulence damping 7^ = 0. The turbulence is assumed to be confined in the 
acceleration region; possible turbulence leakage from the acceleration region is compensated 
by the adopted sustained source of the transverse component of the large-scale turbulence. 
Particles, however, can escape from the region through its boundaries because of a large mean 
free path of the particles outside the region. We, therefore, consider a simplified equation 
for S(k, u, t) 



where the expression for the damping rate of large scale turbulence due to particle accelera- 
tion 7 ap = 9D. The parameter 8 was determined (iteratively) in such a way to preserve con- 
servation of the total energy in the system of the turbulence and the particles with account 
for energetic particle escape from the acceleration region. The standard Crank-Nikolson 
scheme providing second order accuracy in time evolution was applied to integrate Eqns. 
CQH©. The spatial transport term in Eq. (JTJ) was approximated as (time depending) es- 
cape time N/T esc , where T esc = R 2 /4x, R is the characteristic size of the acceleration region. 
Actually, the temporal evolution of \ is very slow since in Eq. (j2j) it is dominated by the 
transverse component of the turbulence T(k, u, t), which decays slower than the longitudinal 
component S(k,u,t). 

The temporal evolution of the turbulence is, thus, solely due to particle acceleration 
effect. Particle acceleration time in the model is longer than turn-over time of large scale 
MHD motions providing u > D. Our simplified model assumes, therefore, that the turbu- 
lence is primary produced in the form of transverse motions with the scale about 2n/ko with 
a gaussian spectrum, which produces the corresponding longitudinal turbulence due to the 



dS(k,u,t) 
dt 



jsrT(k,uj,t) -j ap S(k,uj,t), 



(5) 
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mode coupling (term ^sTT{k,uj,i)) in a system of a finite scale size R (where Rko > 1). 
Furthermore the model accounts only for the evolution of large scale (energy containing) 
motions of k\(p) <C 1, where X(p) is particle mean free path due to scattering by small 
scale (resonant and non-resonant) magnetic field fluctuations; the corresponding diffusion 
coefficient is k = y\(p)/3. We did not consider here the tur bulence cascade to resonant 



small) scales (c.f. iMiller et al.lll996l ; iPetrosian fc Bykovl 120081 ) . Instead, we fixed the "mi 



croscopic" diffusion coefficient n(p) due to small scale fluctuations, provided perhaps by the 
whistler waves, and considered the case of intensive large scale turbulent motions provided 
k(p) < k^ 1 ■ V < u 2 >. The kinetics of particles satisfying this inequality is determined by 
turbulent advection and so does not depend on the details of the "microscopic" diffusion 
coefficient. 

The energy range, where this inequality holds, does depend on the charged particle 
mean free path X(p), which is ultima tely defined by generally unknown level of the small- 



scale turbulence. iBastian et al.l (120071 ) determined the mean free path of the radio emitting 
electrons (a few MeV) to be about 10 7 cm from the characteristic decay time of the radio 
light curves for the case, when the turbulent transport of the particles was independently 
confirmed. This estimate is consistent with the small-scale (tens of cm) turb ulence level 



of ab out ~ 10 7 — 10 5 derived from the decimetric continuum burst analysis (INita et al. 



20051 ). Thus, to be conservative, in the case of electrons our approximation is firmly justified 
at least up to the energy about 1 MeV, where the particle transport is fully driven by the 
large-scale turbulence and as so it does not depend on the actual momentum dependence of 
the mean free path X(p). Then the model accounts for a nonlinear backreaction of acceler- 
ated particles on large scale motions only. Although the assumed presence of the small-scale 
turbulence implies the possibility of the stochastic resonant acceleration along with consid- 
ered here non-resonant acceleration, we do not take into account the resonant acceleration 
explicitly, because the energy density of the small-scale turbulence is much smaller than that 
of large-scale turbulence. As has been discussed, the resonant acceleration can, nevertheless, 
contribute to the injection term Fi(p) in Eq. (TjQ). 

We consider injection of non-relativistic particles of a momentum po, i.e. F\(p) oc 5{p — 
Po). It is convenient to characterize the injection efficiency by the injection energy loading 
parameter 

2 J e(p)F i (p)p 2 dp 
* D(0) p < u 2 > ' W 

where e(p) is the particle kinetic energy expressed via its momentum p. In Figure [T] we 
show the particle distribution function (normalized oc p 2 N) calculated for the model. We 
assumed continuous injection of mono energetic particles (electrons and protons i = e, p) 
with the injection energy loading parameters ( e = 10~ 3 (left panel) and ( e = 0.1 (right 
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panel). 



Although there are apparent differences in particle spectra for different ( e , all our runs 
display clearly soft-hard-soft behavior of the spectra of accelerated particles. The origin 
of this spectral evolution is easy to understand within the proposed model. Initial phase 
of the acceleration occurs in the linear regime (test-particle approximation is still valid 
on this stage), which results in effective particle acceleration by the longitudinal large scale 
turbulent motions and spectral hardening. However, fast particles accumulate a considerable 
fraction of the turbulent energy on this stage and start to exhaust the turbulence, thus, the 
efficiency of the acceleration decreases, which first affects higher energy particles resulting 
in the spectrum softening. 

Another important point, which must be noticed from the figure, is that the slope of 
the spectrum at the late decay phase (red solid curves) depends strongly on the injection 
efficiency ( e . In fact, the spectrum is much steeper in case of strong injection. In practice, 
the spectrum in the right panel is so steep that it is probably undistinguishable against 
background thermal particle distribution. This means that the sequence of the (dash-dotted) 
spectra of accelerated electrons in the right panel will reveal itself as SHH evolution of 
the HXR spectrum. This conclusion is consistent with the fact that the SHH evolution is 
observed in stronger events, where enhanced injection of the charged particles (e.g., protons) 
is likely, a nd with a recent find ing of gradual transitions between SHS and SHH evolution 
fragments (jGrigis & Benzll2007l ). which requires a common acceleration mechanism for both 
SHS and SHH evolution patterns, even though additional spectral hardening in the gamma- 



ray r ange can occur due to relativistic particle trapping in the coronal loops (IKrucker et al. 
20081 ). 



Besides the general SHS evoluti on, we should note that in agreement with previous 
studies of the stochastic acceleration (IMiller et al.lll997l ; iGrigis fc Benal2006l ) the spectra do 
not obey power-laws exactly: break-up and break-down turning points are evident from the 
plots. It should be noted here that since the nonlinear effects were taken into account in the 
model the distribution function calculated for mono-energetic injection will not have any of 
the general properties of the Green function of a linear system. Therefore, one can not build 
the distribution function in the nonlinear case using the superposition principle any longer. 
Nevertheless, the initial stage of the particle acceleration occurs in the linear regime if the 
loading parameter ( e is smaller than unity. Thus, broadly speaking, the general behavior of 
particle spectra evolution as it is illustrated in Figure [I] will hold for other relatively narrow 
(as the first blue curves in Figure [T]) initial particle distributions with a similar loading 
parameter defined by Eq. (jSJ). 



Grigis &; Bend ( 120061 ) demonstrated that the HXR spectra from the thin-target coro- 
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nal sources are the most directly linked to the energy spectra of the accelerated electrons, 
while the properties of the thick-target foot-point sources can essentially be affected by the 
transport effects. Accordingly, we computed the evolution of the thin-target HXR emission 
generated by the evolving ensemble of the accelerated electrons (as in Fig. 1) and then 
derived the evolution of the HXR spectral index at E = 35 keV to c ompa re with the obser- 



vations of the coronal HXR sources reported by iBattaglia fc Benzl ( 120061 ). The theoretical 
dependences of the HXR spectral index on time are shown in Figure [2] by three curves with 
different ratio of the acceleration time to the escape time. The asterisks in the same plot 
show the evolution of the HXR spectral index observed for the coronal source in the Dec, 



04, 2002 event (IBattaglia &: Benzll2006l ). Even though no theoretical curve is the exact fit to 
the data, one can clearly note important similarities between theoretical and observational 
curves including the main SHS behavior and some hardening at the later stage. 

Since the spectral index analysis of the coronal source can in principle be biased by much 
stronger foot-point contribution, a more reliable way of the thin-target HXR analysis could 
be the study of the occulted X-ray flares. However, the thin-target HXR flux is typically weak 
from the occulted coronal sources, so the systematic statistical stud y of the occulted flare s 



reports only the spectral data at around the peak time of the flare (IKrucker fc Linl 120081 ) . 
In some cases, however, it is still possible to derive information on the spectral evolution of 
the occulted flares by integrating the signal during rise, peak, and decay phases respectively. 
An example of the corresponding spectral index evolution in an occulted 06 Sept. 2002 flare 
is shown by three long horizontal dashes in Figured (E. Kontar, private communication). 
The SHS evolution is evident in this instance as well. 

In addition to the mentioned similarities between the theoretical and the observed spec- 
tra, there are also apparent differences. We have to note, however, that the differences 
between the theory and observations are less significant than the difference between the 
spectra observed from different events. Thus, we can ascribe these differences to the varying 
geometry of the source and/or to different regimes of the turbulence generation, cascading, 
damping, and escape, i.e., to those details of the model, which have not been specifically 
addressed within this letter. 

To summarize, we note that taking into account the nonlinearity, which is necessarily 
present in a system where efficient acceleration by strong turbulence occurs, offers a plausi- 
ble way of interpreting both kinds of the characteristic spectrum evolution, SHS and SHH, 
observed from solar flares. A side achievement of the adopted here model of the turbulent 
electron transport is the energy independent escape time from the acceleration region, which 
implies that electrons with different energies leave the acceleration si te simultaneously: the 
property required by measurements of the HXR fine structure timing (jAschwanden 2002 ). A 
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full comprehensive picture of the particle acceleration in flares will require further analysis 
with the shock waves, turbulence cascading, and injection details included, as well as com- 
puting the HXR and gamma-ray spectrum evolution, which we plan to consider elsewhere. 

We are cordially grateful to Marina Battaglia for providing us with the data on the 
spectral index evolution from the coronal sources, and to Eduard Kontar, who obtained and 
provided us with the data on the spectral evolution of the occulted 06 Sept. 2002 flare. The 
work was partly supported by Russian Foundation for Basic Research, grants No. 06-02-16844, 
06-02-16295, 08-02-92228, 09-02-00624 and by NSF grants AST-0607544 and ATM-0707319 
to New Jersey Institute of Technology. We have made use of NASA's Astrophysics Data 
System Abstract Service. 
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Fig. 1. — The temporal evolution of particle distribution function (sequence of p 2 N vs 
p/Po plots, where p/po is the dimensionless particle momentum normalized by the injection 
momentum po) simulated within a flare acceleration region of the scale size R = 147r//co for 
the particle injection energy loading parameters £i = 1CT 3 (left panel) and £j = 0.1 (right 
panel). Particle spectra are shown in 20 logarithmically distributed consequent time frames 
measured in tD(0) starting from 0.01 to 30. For some typical parameters, e.g., R = 2 x 10 9 



cm, 5=300 G, n — 10 9 — 10 11 cm 3 we have va — 2.2 x (10 8 — 10 9 ) c m/s, and the characte ristic 

: 1/D(0) ~ 1-10 s in agreement with HXR (IGrigis fc Benzll2006l ) and 



acceleration time r„ 



radio (IBastian et al.l 120071 ) observations. 
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Fig. 2. — HXR spectral index evolution for theoretically calculated spectra with various 
ratios of the escape time to the acceleration time, T esc /r acc = 5 (solid curve), 1 (dashed 
curve), and 0.2 (dash-do tted curve); and observed from the 04 Dec. 2002 flare, asterisks, 
(IBattaglia &: Benzl 120061 ). and from the occulted 06 Sept. 2002 flare, horizontal dashes, 
(E.Kontar, private communication) . 



